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Abstract. 

In this work we present fully differential cross sections (FDCSs) calculations using post 
and prior version of CDW-EIS theory for helium single ionization by 100 MeV C 6+ 
amu -1 and 3.6 MeV amu -1 Au 24+ and Au 53+ ions. We performed our calculations 
for different momentum transfer and ejected electron energies. The influence of 
internuclear potential on the ejected electron spectra is taken into account in all cases. 
We compare our calculations with absolute experimental measurements. It is shown 
that prior version calculations give better agreement with experiments in almost all 
studied cases. 
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1. Introduction 

The study of electron emission spectra in ion atom collisions has been a field of 
intense activity for years (Stolterfoht et al 1997). For intermediate to high energy 
single ionization there has been considerable theoretical efforts focused in the so-called 
two centre electron emission (TCEE) (Fainstein et al 1991, Pedersen et al 1990). 
Improvement in the description of the ionized electron moving in the presence of both 
residual target and projectile fields after the collision (final state) has been key for the 
correct description of experimental data (Gulyas et al 1995). 

Within distorted wave approximations, it has been shown that, at least for high 
impact energy and multiply charged projectiles, CDW theory of Belkic (1978) used 
together with an appropriate description of the initial bound and final continuum 
electron states, yields best results for doubly differential cross sections (DDCSs) (Gulyas 
and Fainstein 1998, Ciappina et al 2003). However, when the projectile impact velocity 
decreases, the CDW-EIS theory of Crothers & McCann (1983) gives better results, its 
only difference being the choice of the initial state. Moreover CDW-EIS approximation 
is formally free of criticisms regarding the initial state proper normalization, and the 
transition amplitudes have not the divergent behavior that CDW exhibits (Crothers 
1982) (although it has been demonstrated that CDW amplitudes are integrable and its 
DDCSs are well behaved) (Dube and Dewangan 1995). 

The field has experienced a renewed interest as a result of the development of 
the experimental technique known as COLTRIMS (cold-target recoil-ion momentum 
spectroscopy) (Moshammer et al 1994). With COLTRIMS, the projectile's tiny 
scattering angle can be obtained indirectly by measuring the ionized electron and recoil 
ion momenta. Fully differential cross sections (FDCS) for ion impact ionization can be 
measured now and this constitute a challenging ground for existing theories (Foster et 
al 2004). 

The first measurements of the FDCS, for various momentum transfers and ejected- 
electron energies, were reported in 2001 by Schulz et al for single ionization of helium by 
100 MeV amu -1 C 6+ . Theoretical results for this process were made later by Madison et 
al 2002, using several approximation schemes. They obtained reasonable good agreement 
between experiment and theory in the scattering plane for intermediate momentum 
transfer, but the theories used was not able to reproduce the measurements for large 
values of momentum transfer and out of scattering plane. 

Subsequently, experiments with other projectiles and energy ranges have been 
performed. Fischer et al (2003) have reported absolute experimental measurements 
for 2 MeV amu -1 C 6+ single ionization of helium in the scattering plane for various 
momentum transfers and ejected-electron energies. Foster et al (2004) have presented 
3DW-EIS results for the single ionization of helium by 3.6 MeV amu -1 Au 24+ and 
Au 53+ ions. The 3DW-EIS model is a modified version of the CDW-EIS approximation 
and, although the authors found good agreement with 2 MeV amu' 1 C 6+ data, the 
theory did not yield a significant improvement for higher charged ions. Theoretical 
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results calculated using a CDW-EIS model exhibited differences between experiment 
and theory on an absolute scale for emission in the scattering plane, defined by the 
plane containing the initial and final projectile momenta. Their calculations were made 
using a post version of the CDW-EIS theory and an active electron approximation 
with hydrogenic wavefunctions for the initial and final electron states. Indeed, the 
simplest description for the He bound initial state is to assume it has one 'active' and 
one 'passive' electron and that the 'active' electron can be described as moving in the 
effective Coulomb field of the atomic core with an effective charge chosen either: (a) to 
reproduce the ionization energy or (b) so that the continuum wave is orthogonal to the 
initial state. 

A more sophisticated description involves the use of full numerical Hartree-Fock 
wave functions for both initial and final states of the active electron (Gulyas et al 1995, 
Gulyas and Fainstein 1998, Foster et al 2004). Hartree-Fock description, however, does 
not include proper angular correlations in the initial state, and for large perturbations, 
there might be the chance that the projectile interacts with more than one electron in 
a single event. An explicit two-electron description, i.e., a full-blown four-body theory 
for the collision process might be necessary in that case. We have shown that by using 
the prior version of CDW-EIS together with an appropriate Roothan-Hartree-Fock 
description of the initial state and an effective charge coulomb wave function for the 
target - electron continuum, we get similar results to those obtained by using numerical 
Hartree-Fock wave functions (Ciappina et al 2004), for ion impact helium ionization 
DDCSs. 

The aim of this paper is to present post and prior CDW-EIS calculations with 
internuclear interaction between the projectile and the target (N-N interaction) taken 
into account for ion helium single ionization FDCSs at different perturbation regimes. 
Atomic units are used throughout unless otherwise stated. 

2. Theories 

We regard He single ionization as a single electron process and assume that (i) the 
initial state for the 'active' electron is described by a semi-analytical Rothan-Hartree- 
Fock scheme using a 5 parameters wave function (Clemente and Roetti 1974) and (ii) 
in the final state the 'active' target electron moves in the combined Coulomb field of 
the target core with an effective charge Z e ff = 1.6875. The electron-projectile relative 
motion are represented in a CDW-EIS approach, i.e. one eikonal phase in the entrance 
channel and a pure Coulomb distortion in the final one. N-N interaction is treated as a 
pure Coulomb interaction between the projectile with a charge Zp and the true target 
core charge, Zt — 1. 

N-N interaction is taken into account in the transition amplitude a^(p), in the usual 
semi-classical or eikonal approximation, through its multiplication by a phase factor 
(McCarroll and Salin 1978), which for pure coulomb internuclear interaction results in 
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(Crothers and McCann 1983) 



a,: 



Hfy p) = i(pv) 2w a if (p) (1) 

were v = ZpZx/v, v is the velocity of the impinging projectile and p is the impact 
parameter (p • v = 0). (p) (a'^ (p)) is the transition amplitude with (without) 
internuclear interaction. Using two-dimensional Fourier transforms we have for the 
transition amplitude elements, CDW-EIS transition matrix can be written alternatively 
SIS db function of the momentum transfer: 

W) = 7^2 I df) T if (r,) [ dp p 2 - e^') " a if (p) (2) 



(2tt, - 

We solve the integral over impact parameter analytically to obtain: 

W = v h,2i Jf iU J dr,' TM) \n - nT 2 ^ (3) 

The remaining integral in (j3J) is evaluated numerically with an adaptive integration 
scheme. This approximation is valid as long as (i) the projectile suffers very small 
deflections in the collision and (ii) the velocity of the recoil ion remains small compared 
to that of the emitted electron. 

Within CDW-EIS, Transition amplitude can be computed as 



rp+CDW-EIS _ L^-CDW 

in its post version or 



W} 



xt EIS ) (4) 



CDW-EIS _ L-CDW 



[xj UDW W^r 18 ) (5) 

in the prior version, where the initial (final) state distorted wave xt (xj) is an 
approximation to the initial (final) state which satisfies outgoing-wave (+) (incoming- 
wave (— )) conditions. For the initial state the asymptotic form of the Coulomb distortion 
(eikonal phase) is used in the electron-projectile interaction together with a semi 
analytical Rothan-Hartree-Fock description for the initial bound-state wavefunction 
(Clementi and Roetti 1974) 

xt EIS = (2n)- 3/2 exp (iK, ■ R T ) i> ls (v T )S+(r P ) (6) 

where £+(rp) is 

£+(r P ) = exp f-i-^- ln {vr P - v • r P )^ . (7) 

The final state wavefunction is collected into the form (Rosenberg 1973, Garibotti 
and Miraglia 1980, Crothers and McCann 1983) 

x fDW = (27r) -3/2 exp (iK/ . R T ) XtMC~ p {v p ) (8) 

where Cp represents the Coulomb distortion of the ejected electron wave function due 
to the projectile coulomb potential. 

Cp{vp) = N(u P ) iFx (-iis P , 1, -ikprp - ik P • r P ) (9) 
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being vp = the Sommerfeld parameter, and kp is the relative momentum of the 
e-P subsystem. The is the Kummer function and N(vp) is the usual normalization 
factor 

N(v P ) = T(l - iu P ) exp(-7Tz/ P /2) (10) 

being T the gamma function. On the other hand Xt( t t) is the wave function for the 
ejected electron in the field of the target residual ion. 

Xt( t t) = (27r)~ 3//2 exp (iky • rp) N(ut) i-Pi ( — ii*r, 1, — ikprp — ikp • rp) (11) 

being 2/p = j£ and now kp is the relative momentum of the e-T subsystem. We use 
Z T = Z e ff = 1.6875 to model the screened target residual ion as a pure Coulomb 
potential. 

The perturbation potentials Wf in equation (JJJ) and W{ in © are defined by 

(H f - E f ) X J = W fX J (12) 

and 

(H i -E i )xt = W iX t (13) 

where Hf (Hi) are the full electronic final (initial) Hamiltonian (neglecting the total 
center of mass motion) and Ef (Ej) are the total final (initial) energy of the system in 
the cm frame respectively. 

The explicit forms of these operators can be written (Crothers and Dube 1992) 

Wf = -V rT • V rp (14) 



and 



1. 



Wi = -V 2 rp - V rT • V rp (15) 

In the centre of mass frame, the FDCS in energy and ejection angle of the electron, 
and direction of the outgoing projectile is given by (Berakdar et al 1993, Inokuti 1971, 
Bethe 1930) 



where iV e is the number of electrons in the atomic shell, \i is the reduced mass of the 
projectile-target subsystem, Ki (Kf) is the magnitude of the incident particle initial 
(final) momentum. The ejected-electron's energy and momentum are given by E k and 
k respectively. The solid angles dflx and dflk represent the direction of scattering 
of the projectile and the ionized electron, respectively. We use non-orthogonal Jacobi 
coordinates (rp, rp) to outline the collision process. These coordinates are the position 
of the active electron with respect to the projectile (rp) and to the target ion (rp) 
respectively. Also the coordinate Rp is needed, that represents the position of the 
incoming projectile with respect to the center of mass of the subsystem e-T. If we 
neglect terms of order 1/Mp and 1/Mp, where Mp is the mass of the target ion nucleus 
and Mp is the corresponding to the incident heavy ion, we can write Rp = rp — rp. 
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We have replaced the transition matrix in the post and prior schemes (equations 
(jlj) and ©) in the definition of FDCS (fT6*j) and we have applied it to several single 
ionization processes. 

3. Results 

We have performed calculations for different projectiles, spanning a large range of 
perturbation strengths as measured by charge to velocity ratio i] = Zp/v. In figure 
1 we present results for 100 Mev amu -1 C 6+ (Schulz et al 2001) single ionization of 
Helium calculated in prior CDW-EIS, for different values of electron emission energy 
(E e ) and momentum transfer (q = K, — Kj). Calculations are in very good agreement 
with available experimental results. 

In figure 2 we layout results for other theories applied to the same process for 
an intermediate value of electron energy and momentum transfer. We see that prior 
CDW-EIS gives the best results. Even when r? = 0.1 both FBA and post CDW-EIS 
fail to accurately describe the experimental results, although they broadly reproduce 
the angular distribution. 

Figures 3 and 4 show results for 3.6 MeV amu -1 Au 24+ impact ionization of He 
(Fischer et al 2003), calculated in prior and post CDW-EIS. For E e = 4.0 eV, results 
for prior version are in reasonable agreement with experiment, However both theories 
fail to correctly reproduce the strong forward emission peak, which is due to the strong 
projectile electron post collisional interaction (PCI). Trend is similar for E e = 10.0 eV, 
(Figure 4) where we see a better performance in prior version calculations, in particular 
in the prediction of the direct peak position. Note that no renormalization factor is 
included in these calculations. 

In figures 5 and 6 we show prior and post CDW-EIS calculations for 3.6 MeV amu -1 
Au 53+ impact ionization of He. Even when we are stretching the validity range of the 
perturbative treatment a bit too much (rj »s 2, 4.4 for Au 24+ and Au 53+ projectiles 
respectively), angular structure with only one strong peak is correctly predicted in 
prior version while post version of the theory predicts two distinct direct and recoil 
peak. However the position of the peak is not correctly given in prior version, again 
because the theory underestimates the strong PCI between the impinging ion and the 
ejected electron, which shifts the emission towards the forward direction. Both versions 
including N-N interaction fail to yield the correct order of magnitude of experimental 
data. Large projectile charges are likely to induce quite a large polarization in the target. 
Effective charges both for residual target-electron and N-N interactions are probably 
not the same than for lower charged projectiles, and it is indeed very probable that the 
effective charge approach is not a good approximation here. Model potentials taking 
into account polarization effects need to be considered for the target, at least in the exit 
channel, but most probably in both initial and final states. 



7 



4. Conclusions 

We have performed FDCSs calculations for highly charged ion impact ionization of 
Helium. We employed prior and post versions of CDW-EIS theories taken into account 
N-N interaction but otherwise using as simple an approach for electronic wave functions 
as possible. Indeed, use of prior version helps us to avoid the need of more precise wave 
functions for the initial or final electronic state. We found reasonably good agreement 
with experimental data, even for projectile charges for which the system is arguably 
outside the range of validity of a perturbative theory. 

We see that for emission in the collision plane, three body dynamics seems to be 
enough to explain most of the structures observed for low energy emission and low 
projectile charge. For Au 24+ and Au 53+ projectiles the larger emission in the forward 
direction is not well reproduced by the theory but, as said before, those cases are outside 
the range where perturbative treatments are known to be valid. However, if the effect 
of target polarization in the entrance channel and the inclusion of higher orders in the 
exit channel distortions, are taken into account, perturbation based calculations could 
probably be brought closer to experimental results. 
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Figure captions 



Figure 1. FDCS for 100 MeV amu 1 C 6+ single ionization of Helium calculated in 
prior CDW-EIS: solid line; experimental data, (Schulz et al 2001) solid circles, (a) 
E e = 6.5 cV, \q\ = 0.88 a.u. (b) E e = 17.5 eV, \q\ = 1.43 a.u. (c) E e = 37.5 eV, 
\q\ = 2.65 a.u. 



Figure 2. FDCS for 100 MeV amu 1 C 6+ single ionization of Helium calculated 
for E e = 17.5 eV and \q\ = 1.43 a.u. in prior CDW-EIS: solid line; post CDW-EIS: 
dashed line; FBA: dotted line; experimental data: (Schulz et al 2001) solid circles. 



Figure 3. FDCS for 3.6 MeV amu" 1 Au 24+ single ionization of Helium for E e = 4 cV. 
Prior CDW-EIS: solid line; post CDW-EIS: dashed line; experimental data: (Fischer 
et al 2003) solid circles. Note that the angle of electron emission has been changed 
with respect to the other figures and now is measured in a range of —180° to +180°, 
being 0° the direction of the incoming projectile. 



Figure 4. Same as in figure 3 for E e = 10 eV. 



Figure 5. FDCS for 3.6 MeV amu -1 Au 53+ single ionization of Helium for E e =A eV. 
Prior CDW-EIS: solid line; post CDW-EIS: dashed line; experimental data: (Fischer 
et al 2003) solid circles. 



Figure 6. Same as in figure 5 for E e = 10 eV. 
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